Vanishing Meissner effect as a hallmark of in-plane FFLO instability in 
superconductor — ferromagnet layered systems 
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We demonstrate that in a wide class ol multilayered superconductor - lerromagnet structures 
(e.g., S/F, S/F/N and S/F/F') the vanishing Meissner effect signals the appearance of the in- 
plane Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) modulated superconducting phase. In contrast to 
the bulk superconductors the FFLO instability in these systems can emerge at temperatures close 
to the critical one and is effectively controlled by the S layer thickness and the angle between 
magnetization vectors in the F/F' bilayers. The predicted FFLO state reveals through the critical 
temperature oscillations vs the perpendicular magnetic field component. 

PACS numbers: 74.45.+C, 74.78.Fk, 74.62.-c 
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The diamagnetic supercurrent and resulting mag- 
netic field expulsion observed in seminal experiments by 
Walther Meissner [l[ are known to be one of fundamental 
phenomena peculiar to superconducting materials. The 
London theory 0] gives us a famous expression for the 



supercurrent density j 



n s A/mc originating from 



the phase rigidity of the wave function of superconduct- 
ing electrons. Here n s is the density of superconduct- 
ing electrons, m is the electron mass and A is the vec- 
tor potential. Assuming naturally the electron density 
and mass to be positive we always get the j — A rela- 
tion corresponding to a diamagnetic response. Recently 
this observation has been questioned in several theoreti- 
cal works [Jill predicting the sign change in the London 
relation and an unusual paramagnetic response of the 
hybrid superconductor/fcrromagnct (S/F) and supercon- 
ductor/normal metal (S/N) systems. Such anomalous 
Meissner effect has been attributed to the odd-frequency 
spin-triplet superconducting correlations generated due 
to proximity effect 

For S/F systems the inversed sign of the Meissner cur- 
rents is closely related to the oscillatory behavior of the 
Cooper pair wave function inside the ferromagnet 0, [|| ■ 
These oscillations arc known to cause a number of impor- 
tant fingerprints of the S/F proximity effect including lo- 
cal increase in the electronic density of states at the Fermi 
energy 0-12|, 7r— Josephson junction formation [l3. 14| 
and non-monotonic dependencies of the critical temper- 
ature of S/F bilayers on the F layer thickness fl5, 16|. 

The unusual electromagnetic response contribution be- 
comes even stronger for a superconductor placed in con- 
tact with a composite F/F' layer with different mu- 
tual orientations of the magnetic moments. Such sys- 
tems are known to reveal so-called long-range triplet 
superconducting correlations predicted in Refs. H, fl7| . 
The local supercurrent density can be written as j = 
— e 2 (n s — nt)A/mc, where n s {nt) is the density of the 
singlet (triplet) condensate. Different character of the n s 



and nt components decay leads to the change in the sign 
of the local response, i.e., inversion of the Meissner effect. 
During the last two years an important breakthrough in 
the experimental observation of the long-ranged triplet 
proximity effect occurred 18, [l9j]. All this makes very 



timely the study of the magnetic response of the prox- 
imity induced triplet superconductivity. Note that the 
first experimental measurements [20j | of the London pen- 
etration depth in thin S/F bilayers revealed a slightly 
non-monotonic dependence of the penetration depth on 
the F layer thickness, which was in accordance with the 
theoretical analysis [21| . 

In this paper we address the intriguing problem of 
the Meissner response of the S/F systems exhibiting the 
above sign change in the relation between the super- 
current density and vector potential and show that the 
anomalous Meissner effect can cause the in-plane Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) instability 0, H3| 
of the superconducting uniform state. To elucidate our 
main results we start from rather general arguments il- 
lustrating the physical origin of the instability in sys- 
tems with the anomalous Meissner effect. Consider- 
ing the local supercurrent density (j = —SFa/SA = 
—e 2 n s A/mc) as a variational derivative of the free energy 
functional we find the corresponding free energy term: 
Fa = f(e 2 n s A 2 /2mc)dV. The sign change in the cur- 
rent - vector potential relation can be considered as a 
change in the sign of the effective mass. Introducing the 
superconducting order parameter phase ip and writing 
the free energy in the gauge-invariant form 



F A 



2mc 



dV, 



(1) 



where $o is the flux quantum, one can clearly see that 
the negative local effective mass can, in principle, re- 
sult in the instability of the homogeneous state with 
ip = const, A = and appearance of the phase ip modu- 
lation. Namely such a situation is realized at the transi- 
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tion to the non-uniform FFLO state (see the discussion, 
for example, in [IH). As a consequence, the above ex- 
pression describing the linear current response should be 
reconsidered for a new inhomogeneous ground state. 

To illustrate the above general arguments by a concrete 
example of instability we hereafter focus on the consid- 
eration of thin film structures of total thickness much 
smaller than the screening length. This assumption al- 
lows us to consider only the currents flowing in the film 
plane and neglect the change of the vector potential on 
the structure thickness. Introducing the in-plane FFLO 
modulation vector k so that ip = kry we find: 



F A 



S 



e 2 n s 
2mc 



dx , 



(2) 



where the x axis is chosen perpendicular to the film 
plane, S is the sample area in the (yz) plane, An and 
r || are parallel to the film. All the states with A -1 = 
J(e 2 n s /2mc)dx < are clearly unstable and, thus, the 
boundary of the in-plane FFLO instability is given by 
the condition A -1 = of vanishing Meissner effect for 
the in-plane field. Note that the above arguments, be- 
ing applied for the FFLO state itself, clearly show that 
in the modulated state the Meissner response should be 
diamagnetic. Thus, in the systems under consideration 
the paramagnetic Meissner response appears to be im- 
possible. 

We now proceed with the microscopic calculations of 
the FFLO critical temperature and magnetic screening 
length for three particular structures (S/F, S/F/N and 
S/F/F') shown in Fig.Q] Note that for S/F bilayers the 
modulated along the F layer state has been sug gested in 
Ref. [25| but later it has been pointed out |26| that the 
conclusions of Ref. [251 ] are based on the wrong boundary 
conditions assuming the modulation of the order param- 
eter only in the F layer. In contrast with Ref. [25[ in 
our case the same modulation is present both in S and 
F layers. Somewhat similar non- uniform phase has been 
predicted for a ferro mag netic cylinder covered by the su- 
perconducting shell 27|. Interestingly in 3 He films the 
non-uniform superfluid p-wave state may be stimulated 
by the surface scattering of quasiparticles [28| . In our 
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FIG. 1: (Color online) The sketch of the hybrid structures 
under consideration. S layer is placed in contact with (a) 
F film, (b) F/N bilayer and (c) F/F' bilayer with different 
magnetic moment orientations shown by arrows. 



calculations we assume that: (i) the system is in a dirty 
limit; (ii) the exchange field h in the fcrromagnet is much 
larger than the critical temperature T c0 of the isolated S 
layer; (iii) the thickness of the S layer d s is smaller than 
the coherence length £ s = J D s /2irT (D s is the diffu- 
sion constant in a superconductor), so we can neglect 
the variation of the order parameter function A across 
the S layer; (iv) all interfaces are transparent. 

Near the critical temperature the anomalous Green 
function 



fn 

hi 



/12 

/22 



= (f 8 + t t a) io-y. 



satisfies the linearized Usadel equation [30( 

y V 2 / - uj - \ (ha/ + /ha) + A = 



(3) 



(4) 



where A = Aia y is the superconducting gap function, 
Lu n = irT(2n + 1) are the Matsubara frequencies, and 
D is the diffusion constant, which may be different for 
different layers. In the absence of the barriers between 
layers the function / as well as the combination ad x f are 
continuous at each interface (a is the Drude conductivity 
of the corresponding layer) . Wc assume Fermi velocities 
in all layers to be equal, so that the ratio between con- 
ductivities of different layers is the same as the ratio be- 
tween the corresponding diffusion constants. The critical 
temperature T c of the system is determined by the com- 
ponent fi 2 of the Green function in the superconductor 
in accordance with the self-consistency equation 



Aln£+ £ 



T c0 



A 



2n + H 



- 7iT c /£ = 0, (5) 



where T c0 is the critical temperature of the isolated su- 
perconducting film. 

In the limit of weak screening the Meissner response 
averaged over the structure thickness do takes the form 



A = 



Ado ec$nd n ^„./ 



I6tt 3 T c 



dx 



(6) 



This expression clearly shows that the triplet compo- 
nent provides the negative contribution to A~ 2 . To 
describe the FFLO state we assume the gap A(r y ) = 
Aoexp(ikr||) and the anomalous Green function / = 
tp(x) exp ukriij to be spatially modulated. 

We start from the simplest case of a bilayer (see 
Fig. (Ua)), which consists of a thin S film and F layer 
of the thickness df <C The exchange field h in the 
F layer is uniform and directed along the z-axis, so that 
fu = f-22 = 0. Substituting the modulated Green func- 
tion into the Usadel equation and solving it under the 
assumption, that the function / weakly varies across the 
S layer, we obtain the components fi% in the S(F) layer 



J 12 — 
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u n + T s 1 (k) ' 



fl2 



s cosh (q k (x 
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df)) 



cosh (q k d f ) 



(7) 



3 




0.09 0.1 0.11 0.12 0.13 0.14 u o 0.2 0.4 0.6 0.8 1 

(WW % 



FIG. 2: (Color online) Magnetic screening parameter X~ 2 
(blue solid lines) in the uniform superconducting state and the 
optimal FFLO modulation vector ko (red dashed curves) vs 
the S layer thickness d a for the S /F bilayer. Blue dash-dotted 
curves drawn by hand illustrate the behavior of the magnetic 
screening parameter in the FFLO regime while green dotted 
curves correspond to the A -2 behavior calculated for the un- 
stable uniform state. We take £ s o = \J D 3 /4:nT c o — 0.1£/ and 
(1) d f = 0.75&, (2) d f = 1.0O, (3) d f = 1.2C/, (4) d f = 
2.0f/. Also we denote X~f f = A~ 2 (T c (0)ec$ do/27ra s d s A 2 ). 

where q k = yjq 2 + k 2 , q = (1 + i) /£/, and £/ = yjDfjh 
is the coherence length in the ferromagnet. The complex 
pair-breaking parameter 

T-\k) = ^k 2 + ^ a -L qki& nh{q k d f ) (8) 

determines the critical temperature T c (k) of the S film: 

where Vl 7 is the Digamma function. Note that these re- 
sults can be obtained by replacing us n — > cu n + D s ^k 2 /2 
in the Usadel equation for the uniform state. 

The effective magnetic screening length in the uniform 
state can be expressed through the derivative of the above 
expression for T c at k = 0: 

a d s a s A 2 D s \ f (1 All 0T C 

= -2ecWc(0) [^M^U J/J W * 

(10) 

Here v = t s 1 (0) / 2ttT c (0) and "Ji is the trigamma func- 
tion. Calculating the derivative of the Eq. (j9]) we find the 
result obtained in [2l| . The condition of the stability of 
the uniform superconducting state, dT c /dk 2 {k = 0) < 0, 
imposes a diamagnetic character of the Meissner response 
for the magnetic field parallel to the plane of the layers. 

For df ~ £/ the contribution from the F layer to 
the Meissner response coefficient A -2 can become neg- 
ative. For S/F bilayers with a large difference in the 
diffusion constants (Df/D s 3> h/T c o) the screening pa- 
rameter A -2 can even vanish at some critical thickness 



FIG. 3: The dependencies of the critical temperature T c vs 
the modulation vector k f or different thicknesses d s of the 
S film. We take £ s0 = V 'D s / '4ttT c0 = 0.1£/, d f = 1.2^ 
and the following set of values (d a /^f)(a s /af): 0.13, 0.125, 
0.12, 0.1165, 0.1148, 0.114, 0.113, 0.1125. The increasing d s 
thickness corresponds to the increasing T c maximum. 

d sc ~ (a f I '<J s )£f ■ At the critical thickness d s = d sc the 
derivative dT c /d(k 2 ) \ k=0 turns to zero and d s < d sc the 
superconducting transition occurs not to the uniform but 
to the modulated FFLO state with the modulation vector 
ko 0. The typical dependencies A~ 2 (d s ) arc shown by 
blue solid curves in Fig. [5] The dependencies k (d s ) for 
different df are shown by red dashed curves in Fig.O The 
corresponding dependencies T c (k) are shown in Fig. [3] ft 
is interesting that the FFLO state can survive even for 
parameter range corresponding to a complete suppression 
of the uniform BCS state for all temperatures. 

Discussing the physical reason of the FFLO phase 
emerging in S/F bilayer we should note that the prox- 
imity effect with a ferromagnet plays a role of the pair- 
breaking effect at the S/F interface. The FFLO-likc 
modulation of the pairing wave function weakens such 
pair-breaking, but at the same time this modulation sup- 
presses partially the critical temperature of the S layer. 
The efficiency of the first (second) mechanism is gov- 
erned by the diffusion coefficient Df(D s ) and finally for 
Df/D s 3> h/T C Q the critical temperature of the FFLO 
state may exceed that of the uniform one. 

Our analysis reveals a direct relation between the van- 
ishing Meissner effect and the FFLO phase formation. In 
Fig. [4] wc show the distribution of the triplet and singlet 
components over the bilayer thickness in the BCS state 
close to the threshold of the FFLO instability. One can 
see that the triplet component providing the anomalous 
contribution to the Meissner effect strongly exceeds the 
singlet one at the free surface of the F layer. This cir- 
cumstance gives a hint how to stabilize the FFLO phase: 
one should add the normal metal (N) layer on the top of 
the ferromagnetic layer (see Fig. 1(b)). Moreover such 
modified system allows to overcome the strong damping 
of T c in the FFLO state of the S/F bilayer and get the 
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FFLO state for temperatures close to T c q. Details of cal- 
culations can be found in Supplemental Material . 

The appearance of the FFLO state can be effectively 
controlled provided we consider S/F/F' structures (see 
Fig. 1(c)) with a certain angle 8 between the magnetiza- 
tion vectors in the F and F' layers. Such systems are 
recently discussed as possible candidates for spin valve 
devices 3ll - 33[ . For non-collinear magnetic moments the 
triplet component of the anomalous Green function gen- 
erated in the F film becomes long-range in the F' layer 
and decays at a distance of the order of 3> £/ (where 
£ n = y/Df /4irT c0 ) while the singlet component is fully 
damped at a distance ~ = \JDf /h from the F/F' 
interface. As a result, if the thickness df of the F' layer 
strongly exceeds then the corresponding contribution 
into the screening parameter A~ 2 is always negative and 
can become comparable with the one from the S film. 
In the simplest case of small (large) thickness of the 
F(F') layer, i.e. df -C £/ and df — > oo, the long-ranged 
triplet component //^ in the F' layer is proportional to 
(df /£/) 2 sin 8. For large diffusion constant Df the ratio 
between the negative contribution coming from the F' 
layer and positive S layer contribution to the screening 
parameter A~ 2 can become of the order of unity for 



sin 



> 



D s d s 



(11) 



Varying the angle 8 one can trigger the transition from 
the uniform state, realized for 8 close to zero and tt, to 
the FFLO state, which is favorable for 9 close to ir/2. 
Thus, the formation of the FFLO phase should affect the 
angular dependence of the critical temperature in S/F/F' 
spin valves devices. 

Experimentally the appearance of the FFLO state can 
be identified by the observation of the critical tempera- 
ture oscillations vs mag netic field H perpendicular to the 
plane of the layers [34[. For simplicity we consider here 
only the case of a S/F bilaycr. Choosing an appropriate 




FIG. 4: (Color online) The spatial profile of the singlet (blue 
dashed curve) and triplet (red solid curve) components of the 
anomalous Green function in the S/F bilayer at T — T c (0) 
and u„ = 7rT. We take £ s o = 0.1£/, df = L2£/ and 
(d s Ji f ) (a s /a f ) = 0.13. 




FIG. 5: (Color online) The phase diagram of the S/F bilayer 
in the FFLO regime (red curve). Dashed curves correspond to 
dependencies T C (H) for different n. We put here £ s o = 0.1£/, 
df — 1.2£/ and d s = 0.12(o"//a" s )C/- For comparison in the 
inset we show the H-T phase diagram for d s = 0.14(<t//<t s )£/ 
corresponding to the uniform superconducting state. We de- 



note H = $o /47r£y 



vector potential A(rn) in the plane of the layers we get 
the Usadcl equation for the component fi2 in the form 



D 



0. 



9 -„-^A(r||) 



The solution of the Eq. (TT 
Xn (r||) ip(x), where Xn (ry 



fi2-(u„ + ih)fi2+A = 0. 

(12) 

takes the form: /12 — 
is an eigenfunction of the 

Hamiltonian H = — [<9 r|| — 27rA(r||)/$o] 2 - The critical 
temperature corresponding to the n— th Landau level is 
defined by Eq. © with k 2 -t 2irH(2n + l)/$ . The 
competition between levels with different n results in a 
peculiar dependence T C (H) shown in Fig. 4. 

In conclusion, we predict that vanishing Mcissner ef- 
fect in thin-film multilayered S/F systems should result 
in the in-plane FFLO instability which is particularly 
important for designing the 7r— junction or spin valve de- 
vices. Interestingly that in contrast with the original 
FFLO phase [H,^ which emerges at relatively low tem- 
perature, the modulated phase in S/F/N or S/F/F' het- 
crostructures may appear near the critical temperature of 
S layers. The appearance of the FFLO phase besides the 
anomalous behavior of the screening, should also result 
in the oscillatory-like temperature dependence of the per- 
pendicular critical field. For the S layer thickness exceed- 
ing the screening length the FFLO instability can not, of 
course, expand into the bulk superconductor. In this 
case the surface instability can reveal through the forma- 
tion of the vortex sheet parallel (and positioned close) to 
the S/F interface suppressing, thus, the anomalous part 
of the Meissner response. The linear vorticity density 
should be proportional to the FFLO modulation vector 
ko. Note in conclusion that similar instabilities could also 
appear for the superconductors with anisotropic pairing 



5 



where the paramagnetic currents are caused by the sur- 
face - induced Andreev bound states [35, 
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